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ABSTRACT 

We present the detailed optical to far-infrared observations of SST J1604+4304, 
an ULIRG at z = 1.135. Analyzing the stellar absorption lines, namely, the Call H & 
K and Balmer H lines in the optical spectrum, we derive the upper limits of an age for 
the stellar population. Given this constraint, the minimum x 2 method is used to fit 
the stellar population models to the observed SED from 0.44 to 5.8u\m. We find the 
following properties. The stellar population has an age 40 - 200 Myr with a metallicity 
2.5 Zq. The starlight is reddened by E(B — V) = 0.8. The reddening is caused by 
the foreground dust screen, indicating that dust is depleted in the starburst site and 
the starburst site is surrounded by a dust shell. The infrared (8-1000u.m) luminosity 
is Li r = 1.78 ± 0.63 x 10 12 Lq. This is two times greater than that expected from 
the observed starlight, suggesting either that 1/2 of the starburst site is completely 
obscured at UV-optical wavelengths, or that 1/2 of Li r comes from AGN emission. The 
inferred dust mass is 2. Oil. Ox 10 8 Mq. This is sufficient to form a shell surrounding the 
galaxy with an optical depth E(B — V) = 0.8. From our best stellar population model - 
an instantaneous starburst with an age 40 Myr, we infer the rate of 19 supernovae(SNe) 
per year. Simply analytical models imply that 2.5 Zq in stars was reached when the 
gas mass reduced to 30% of the galaxy mass. The gas metallcity is 4.8Zq at this 
point. The gas-to-dust mass ratio is then 120 ± 73. The inferred dust production 
rate is 0.24 ± 0.12Mq per SN. If 1/2 of Li r comes from AGN emission, the rate is 
0.48 ± O.24M per SN. We discuss the evolutionary link of SST J1604+4304 to other 
galaxy populations in terms of the stellar masses and the galactic winds, including 
optically selected low-luminosity Lyman ev-emitters and submillimeter selected high- 
luminosity galaxies. 

Key words: cosmology: observations — dust, extinction — galaxies: individual: SST 
J1604+4304 — galaxies: evolution — galaxies: stellar content — galaxies: starburst 
— galaxies: high-redshift — infrared: galaxies 



1 INTRODUCTION 

IRAS (Infrared Astronomical Satellite) observations discov- 
ered numerous infrared galaxies. At bolometric luminosities 
> 10 11 I/q, infrared galaxies are the dominant population 
in the local Universe, being more nume rous than optically 
selected starbursts, AGNs, and quasars (|Sanders fc Mirabell 



1996). By luminosity, infrared galaxies are classified into lu- 
minous infrared galaxies(LIRGs) with Li r > 10 11 I/q, ul- 
traluminous infrared galaxies(ULIRGs) with Li T > 10 12 Z/q, 
and hyperluminous infrared galaxies(HyLIRGs) with Li r > 
lO 13 Z/0. The ratio of infrared to visi ble luminosity, Lir/Ls 
increases with L ir l|Soifer et al.ll 19891 ). Although there is ev- 
idence that an optically buried AGN may exist in LIRGs, 



E-mail: kkawara@ioa.s.u-tokyo. ac.jp 



i is ongoing dArmus et al.l 
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optical spectroscopic data support the superwind model in 
which the kinetic energy provided by supernovae (SNe) and 
winds from massive stars in starburst site drives a large-scale 
outflow. 

The Infrared Space Observatory (ISO) source counts in 
the mid- and far-infrared extended our knowledge of infrared 
galaxies up to ; ~ 1 ( Taniguc hi et al 
1997:Kawara et al. 
1999 Serj eant et al.l 



1998; 



2000; 



Puget et al.l 



1997: Oliver ct al 



19991; lElbaz et al 



Sato et al ] 120031 : [k awaraetal 



20041 ). and revealed the number of ULIRGs rapidly grew with 



increasing redshift(Gcnzel & Cesarskv 2000; Charv & Elbaz 
l200ll ; lHeraudeau et al.ll2004l ; IOvabu et al.ll2005l 1. 

The next large advance was made by the Spitzer Space 
Telescope with the sens itivity to probe in frared galaxies at 
z = 2 — 4 (se e revie w in lSoifer et al.ll2~008l reference therein). 
iReddv et all (2006), using deep Spitzer 24um observations, 
examined star formation and extinction in optically selected 
galaxies, near-infrared-selected galaxies, and submillimeter- 
selected galaxies (SMGs) at z ~ 2. Spitzer provided a pow- 
erful probe to measure redshifts of SMGs, the most lumi- 
nous i nfrared galaxy population as HyLIRGs dPope et alj 
l2008al ). Based on Spitzer mid-infrared data, the diagnos- 
tic diagram to separate starburst- and A GN-dominated in- 
frared galaxies has been developed(e. g., IStern et"afl 120051; 
iPope et al.ll2008bl ; |PoTletta et al.ll2008h . and it is suggested 
that Li r in SMGs as wel l as ULIRGs at z ~ 2 is dom- 
inated by star formation l|Pope et al l l2008ah . while many 
SM Gs contain an AG N as evidenced by their X-ray proper- 
ties. iDev et all l|2008h uncovered a new population of dust- 
obscured galaxies (DOGs). These galaxies have Li r compa- 
rable to ULIRGs. However, their Lir/Ls is greater than 
local ULIRGs. DOGs and SMGs as well as high-redshift 
(z ~ 1) ULIRGs are generally considered to be young and 
massive galaxies, although the nature of the stellar popula- 
tions is far from being understood. To understand the evo- 
lutionary link of these infrared galaxies to optically selected 
and near- infrared-selected galaxies, it is crucial to observe 
stellar absorption lines which indicate the age of the stellar 
populations. 

During the course of searching for distant objects in a 
cluster field, SST J1604+4304 drew our attention with its 
red optical-3.6p.rn color. Follow-up study revealed this object 
is an ULIRG at z ~ 1.135 and its energy source is young 
stars with an age < 200 Myr. 

We present the data consisting of optical spectroscopy 
and photometry from the optical to far-infrared in section 2. 
We derive an age and metallicity of the stellar population, 
and show a foreground dust screen is plausible in Section 
3. In Section 4, mid- and far-infrared emission by dust and 
energetics are discussed. In Section 5, we discuss the inferred 
size of a dust shell surrounding the galaxy, the SN rate and 
metallicity expected from the broad-band SED analysis, and 
the dust production rate per SN. The evolutionary link of 
infrared galaxies to optically selected galaxies is discussed 
based on the galactic wind models of elliptical galaxies. 

We adopt H = 70 km s _1 Mpc" 1 , Q. m = 0.3, and 
SIa = 1 — fim throughout this paper. In this cosmology, 
the luminosity distance to the object is 7520 Mpc. The flux 
density F v is used in units of /xJy, which is converted to AB 
magnitudes through the relation m(AB) — —2.5log(F v ) + 
23.9. 




Figure 1. Images centered at SST J1604+4304, where the size is 
70" X 70"for 78i4, 3.6pm, 8.0pm, and 24pm, and 240" X 240"for 
70pm, and 160pm. For all plots, north is at the top and east is 
to the left. We present 32" radius circles which are the size of the 
aperture used for MIPS 160 pm photometry. Two bars in the Isn 
image are used to indicate SST J1604+4304. 



2 OBSERVATIONS AND DATA ANALYSIS 

This work is based on the archival data taken by the Spitzer 
Space Telescope, the Hubble Space Telescope (HST), and the 
Subaru Telescope and new photometric observations per- 
formed on the United Kingdom Infrared Telescope(UKIRT) 
and the telescope. In addition, the Gemini 

Telescope North was used to do optical spectroscopy. 



2.1 Target 

SST J1604+4304 has been found at 

R.A. 16 h 04 m 25!538 
deck +43°04'26'.'55, 

in the J2000 system in a field observed with the Spitzer 
IRAC instrument. The optical counterpart in the HST data 
is located within the 0'.'4 offset accuracy requirements of 

1 This telescope is dedicated to the Multicolor Active Galactic 
Nuclei Monitoring (MAGNUM) p roject and is installed at the 
Haleakala Observatories in Hawaii jYoshiill2002T) . 
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Figure 2. 10" X 10"images centered at SST J1604+4304. For all 
plots, north is at the top and east is to the left. We present 3'.' 7 
radius circles which are the size of the aperture used for IRAC 
3.6-8.0um photometry. 



Spitzer (| Werner et al.ll2004l ). Images of SST J1604+4304 
are given in Figures Q] and [2] The broad-band fluxes are 
given in Table [T] 

SST J1604+4304 is only 32" away from the center of 
the massive g alaxy cluster CL 1604+4304 at z = 0.90 



dGal fc Lubinl 2004 ). Based on a weak-lensing analysis, 



Margonincr ct al. riooi) find the mass contained within pro- 



jected radius R is (3.69 ± 1.47) [R/{500kpc)] x 10 14 M Q , cor- 
responding to an inferred velocity dispersion 1004 ± 199 
km s" 1 . Thus, using a singular isothermal sphere, we ob- 
tain the magnification m = 1.17 + 0.09 for SST J 1604+4304 
l|Schneider.P.. Ehlers.J.. fc Falcdfl992l ). 



2.2 Spitzer data 

All IRAC and MIPS data sets were retrieved from the 
Spitzer public archive. The IRAC data consist of four broad- 
band images at 3.6, 4.5, 5.8, and 8.0 um. The field-of-view 
is nearly 5f 2 x 5^2 with a linear sc ale of approx i matel y l"22 
pixel -1 . The reader is referred to Fa zio et al.l {2004) for a 
complete report on the in-flight performance of IRAC. In the 
pipeline processing, ten BCD FITS images were co-added 
into a single mosaic image per band with a technique simi- 
lar to 'drizzle'. The total effective exposure time is 1936 sec 
pixel - per band. Photometry was performed using a 3 pixel 



(3'.'7) radius aperture. As can be seen in the 3.6u.m image in 
Figure [1] there are two sources which can contribute to the 
flux within the photometric aperture; a faint source is at 6" 5 
to the southwest and a bright source 10"to the west. Using 
the point-spread function which is available on the IRAC 
website, their contributions were estimated to be 3.1 - 4.3 
% of the flux of SST J1604+4304 depending on the pho- 
tometric band. The total fluxes after applying the aperture 
correction are given in Table [T] The quoted errors include 
5% absolute calibration uncertainty that SSCfl recommends 
for general observers. Note that the color correction is very 
small and ignored. 

The MIPS data were taken at 24, 70, and 160 um in 
2004 March. In the pipeline processing, the mapped BCD 
data were rebinned and coadded to a single image with a 
linear pixel scale of 2'.'45, 4'.'0, and 4'.'0 pixel -1 for 24, 70, 
and 140 urn, respectively. The total integration times are 
93, 53, and 84 sec pixel -1 . In the following analysis, we use 
the following post-BCD products, namely, mosaic images for 
24um and filtered mosaic images for 70 and 160u.m. 24um 
photometry was performed using a 7" radius aperture. The 
contributions, to the flux within the photometric aperture, 
from the two sources (see the 24^im image in Figure [TJ, lo- 
cated approximately at 15" to the west-northwest and 13"to 
the southwest, were estimated using the point-spread func- 
tion on the MIPS website. Their contributions are 9% of the 
flux within the aperture. Then, the color correction for 500K 
blackbody and the aperture correction were applied. 

MIPS 70 and 160 urn photometry was performed using 
16"and 32"radius apertures, respectively. As seen in Figure 
CD SST J1604+4304 is faint in the far-infrared. We mea- 
sured sky fluxes with the same aperture at 10 or 11 po- 
sitions along the annulus with radii of 36"- 72"at 7Q\im 
and 96 - 150" at 160u.m. The flux averaged over these sky 
fluxes was subtracted from the flux measured within the ob- 
ject aperture centered at SST J1604+4304. The aperture 
correction and the color correction for 30K were applied. 
The detections are 2.5cr and 5. Oct at 70 and 160u.m, respec- 
tively. As seen in the 24u.m image in Figure [T] there are 
several 24u.m sources within the 160u.m photometry aper- 
ture. Judging from their [3. 6 urn] -[24 urn] color and 70 and 
160 \im brightness at their respective positions, the source 
at 13'.'5 to the west-northwest, the only source as red as 
SST J1604+4304, may contribute 1/3 of the flux within the 
aperture at most. The 160 urn flux is 38.4 ± 7.65 mjy if 
the contributions from this source can be ignored, while the 
flux is 25.6 ± 7.65 mjy if this source contributes 1/3 of the 
flux within the aperture. The real flux of SST J1604+4304 
should be between them, and thus the flux is 32 ± 14 mjy. 
The aperture correction and the color correction for 30K 
blackbody were applied. Note that the 24u.m error in Table 
[1] includes 4% absolute calibration uncertaintjjf). 



2.3 Optical imaging data 

The V606, ^814, and Zsso image data were retrieved from the 
HST public archive, where V606, Int, and Zsso denote the 



2 See Infrared Array Camera Data Handbook 

3 Multiband Imaging Photometer for Spitzer (MIPS) Data Hand- 
book 
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F606W, F8UW, and F850LP, respectively. The Advanced 
Camera for Surveys (ACS) Wide Field Camera (WFC) was 
used, which covers a field of 3l4 x 3l4 with a scale of 0'.'05 
pixel -1 . Four image data were combined into a single co- 
added image per band. The total exposure times are 4840 
sec pixel -1 for Veo6 and Jgi4 and 6000 sec pixel -1 for Zgso- 
The Isu image in Figure [1 shows that SST J1604+4304 
has irregular morphology extending l'.'6 in the NW direc- 
tion with a 0'.'6 width. To determine the total flux of this 
object, /8i4 photometry was performed using 0'.'8, l'.'O, l'/2, 
and 270 radius apertures. The 7si4 flux grows from 0'.'8 to 
l'.'2, while it barely increase from l'.'2 to 2'.'0. Note that con- 
tributions from the two nearest sources were subtracted from 
the flux measured with the 2'.'0 radius aperture. Thus, the 
flux measured with the l'.'2 radius aperture represents the 
total flux of SST J1604+4304. The ratio of the flux with 
a / .'8 aperture to a l'.'2 aperture is 1.24 ± 0.05. Vsoe and 
Zsso photometry was performed using a 0'.'8 aperture. The 
resultant fluxes were converted to the total fluxes by multi- 
plying 1.24 ± 0.05, because we did not measure significant 
changes in optical colors across the target. The total fluxes 
are given in Table [1] The errors in Table Q] include the error 
of aperture conversion from 0'.'8 - l'.'2 and the correction for 
correlated noise in drizzling. 

Subaru B,V,R c ,I c ,z data were retrieved from the 
SMOKA science achieve. The data were taken with Suprime- 
Cam (SUP) imager which has a field-of-view of 34' x 27'with 
a scale of 0'.'2 pixel -1 . The total exposure time is 2880, 
2520, 3600, 1680, 3300 sec pixel -1 after co-addi ng 4 to 11 
imag e frames, using the standard SUP script (jYagi et al.l 
2002). Flux sc aling was performed o bserving 11 standard 
stars given by iMaiewski et ail (I1994T ). I c photometry were 
performed using three apertures with a radius of l'.'O, l'.'2, 
and 2'.'0, and it was confirmed the Subaru I c brightness dis- 
tribution is identical to that measured at the HST 7gi4 pho- 
tometry band. Thus, we use the l'.'2 radius aperture centered 
on SST J1604+4304 to obtain the fluxes given in Table □ 
The errors in Table [T] include 5% photometric uncertainty. 



Table 1. Multiwavclcngth photometry of SST J1604+4304. 



T7V 1 J __ CI 

b liter 


Total flux 


r> b 


Instrument 


UT 


band 


uJv 


arcsec 




vv/mm 


Veo6 


0.224 ± 0.062 


1.2 


HST ACS 


03/08 


^814 


1.05 ±0.08 


1.2 


HST ACS 


03/08 


•^850 


1.53 ±0.13 


1.2 


HST ACS 


07/02 


3.6|jm 


60.35 ±3.16 


3.7 


Spitzer IRAC 


04/03 


4.5|am 


53.76 ±2.79 


3.7 


Spitzer IRAC 


04/03 




£^11 A- A 07 


7 




rtA /f\'i 
U4t/ Uo 


8.0um 


65.68 ±5.36 


3.7 


Spitzer IRAC 


04/03 


24 lim 


319 ± 18.4 


7.0 


Spitzer MIPS 


04/03 


70 nm 


2150 ± 870 


16 


Spitzer MIPS 


04/03 


160 nm 


32000 ± 14000 


32 


Spitzer MIPS 


04/03 


B 


0.0803 ± 0.013 


1.2 


Subaru SUP 


00/06 


V 


0.111 ±0.015 


1.2 


Subaru SUP 


05/05 


R c 


0.218 ±0.018 


1.2 


Subaru SUP 


01/04,05 


Ic 


0.689 ± 0.042 


1.2 


Subaru SUP 


01/04 


Z 


2.29 ±0.21 


1.2 


Subaru SUP 


01/06 


11 


5.40 ± 1.19 


1.2 


UKIRT UIST 


06/02 


K 


20.06 ±3.35 


1.2 


MAGNUM MIP 


06/08 


20cm 


< 0.66 mjy c 


2.7 


VLA e 


1995 


0.2-2 keV 


< 2.1 d 


15 


XMM S 


02/02 


2-10 keV 


< 15 d 


15 


XMM? 


02/02 



a V & o6,hi4.,^Z S 5o denote F606W,F814W,& F850LP, 
respectively. 

b The radius of the photometric apertures which were used to 
obtain the total fluxes. 
c To 3 a. 

d To 3 a in units of 10 -15 erg cm -2 s - 1 . 

e From the FI RST survey jWhite et al.lll997l) . 

f Observed bv lLubin. Mulchaev. fc Postmanl J2004h . The FLIX 

online server was used to derive the upper limit. 



2.4 H/K imaging observations 

The ff-band image was taken with the UIST (1-5 micron im- 
ager spectrometer) on UKIRT in February 2006. UIST was 
used in the imaging mode with a field-of-view of 2'x 2' and 
a scale of 0'.'12 pixel -1 . The total exposure time is 1800 sec 
pixel -1 after co-adding individual image data with a integra- 
tion time of 60 sec. The sky condition was photometric with 
0'.'6 seeing. The K-band image was taken with the multicolor 
imagin g photometer (MIP) on the 2 m MAGNU M tele- 
scope i|Kobavashi et all [l99SllMinezaki et all \2004 ) in Au- 
gust 2006. MIP has a field-of-view of 1.'5 x 1.'5 with a linear 
scale of 0'.'346 pixel -1 . The total integration time is 2340 sec 
pixel - after co-adding individual image data. The weather 
condition was photometric with l"seeing. Flux calibration 
was performed using 2MASS 16042631±4303413 (H=14.36, 
K= 14.39) observed simultaneously with SST J1604±4304. 
Photometric uncertainty is 5% and 9% for H and K, respec- 
tively. 



2.5 Optical spectroscopy 

Optical spectroscopy was carried out using the Gemini 
Multi-Object Spectrograph (CMOS) on the Gemini-North 
telescope in 2007 August. The R150_G5306 grating was used 
along with the OG151_G036 filter, covering 5,000 - 10,000 
A simultaneously with 6.96 A/pixel after binned by 4 in the 
dispersion direction. The different wavelength centers (8,200 
and 8,250 A) were set to fill the CCD gaps. A l'.'O width slit 
was set so that the dispersion direction is perpendicular to 
the long axis of the target. Seven 1800 sec exposures were 
taken in the Nod & Shuffle mode. Data were reduced in a 
standard manner using the IRAF GEMINI. GMOS package. 
The sky subtraction was made with the GNSSKYSUB rou- 
tine. Wavelengths were scaled with CuAr spectra and the 
relative photometric calibration was performed using the 
standard star (HZ44). The results are shown in the bottom 
panel of Figure [3] The spectrum which is median-smoothed 
with a width of 40 pixies (i.e., ~ 280 A), is compared with 
the photometry from the images in the top panel of Figure 
[3] The two measurements agree with each other. 
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Figure 4. The D(CaII) index as a function of ages of stellar 
populations. The index is defined as D(CaII) = 2 X EW(CaII 
K)/EW(H c +CaII H + H e ). The metallicities are 0.005 Z & (long 
dashes), 0.02 Zq (dash dot dot dot), 0.2 Zq (dash dot), 0.4 Zq 
(dashed), 1.0 Zq (dotted), and 2.5 Zq (solid). The three hori- 
zontal bar indicate the observed upper boundary of the D(CaII) 
index at 1, 2, & 3 a upper limits. 



Figure 3. The top panel shows UV to far-infrared SED for 
SST J1604+4304. The broad-band data arc shown along with 
the median-smoothed optical spectrum (gray) taken with Gemini- 
North GMOS. The bottom panel shows the optical spectrum (top) 
together with the synthetic spectrum (bottom) at an age of 400 
Myr, reddened by E(B-V) = 0.6 



3 STELLAR POPULATIONS AND DUST 
EXTINCTION 

3.1 Spectroscopic diagnosis 

The redshift of SST J1604+4304 is 1.135, as identified in 
Figure [3] The synthetic spectrum is also shown for compari- 
son. Th is spectrum was made usi ng the instantaneous-burst 
model (|Bruzual fc Charlo j I2003T ) at an age 400 Myr, red- 
dened by E(B-V) = 0.6. The UV-optical spectrum is char- 
acterized by stellar absorption lines, which are H E at 3970 
A, H<; at 3889 A, and probably higher order lines in the HI 
Balmer series as well as possible presence of the Call HK 
absorption lines at 3968 and 3934 A. The emission line of 
[Oil] at 3727 A is detected. No lines indicating the presence 
of an AGN, such as MgIIA2798 and [NeV]A3346,3426, are 
detected. The 3 a upperlimit of the Mgll EW (equivalent 
width) is 22 A, assuming a FWHM 6000 km s" 1 for the line 
width. This upper limit is small compared to 52 A which 
is the Mgll EW of the co mposite spectrum obtained from 
184 quasi-stellar object bv lTelfer et al.l (|2002l '). There are no 
signs indicating AGN emission at other wavelengths includ- 
ing X-ray and radio as shown Table[T] Thus, we assume that 
the UV-optical spectrum is dominated by starlight. 

The amplitude of the 4000 A discontinuity is known as 
a useful tool to measure the age of stellar populations(eg., 
lBruzuaj[l983l ;l Matsuoka et al.|[2008l ). However, it is not sim- 



ple to correct for the reddening for heavily obscured objects 
like SST J1604+4304. Alternatively, the equivalent widths 
of Call K a nd Balmer H lines are used to sca le ages of 
star clusters (|Rabinlll982l ; ISantos fc Piattill2009h . Here We 
use the index which is defined as D(CaII) = 2 x EW(CaII 
K)/[EW(H c )+EW(CaII H + H e )], where EW(CaII K) is 
the equivalent width for the Call K line, EW(H^) for the 
H c line, and EW(CaII H + H E ) for the blends of the Call 
H & H e lines. The equivalent width is obtained by integrat- 
ing across the line as EW = f x ?(F\ — F\)d\/F c , where 

F° = J™F£d\/{\2 - Al). (A1,A2) in A is (3871, 3909), 
(3926, 3940) and (3951, 3990) for H c , Call K, and Call H 
+ H e , respectively. F\ is the observed flux density and F£ is 
the continuum across the line, which is derived by fitting a 
linear function to assumed continuum points on either side 
of the line profile. Thus, we have D(CaII) = 0.03 ± 0.11 
for SST J1604+4304. where EW(CaII K) = 0.4 ± 1.4 A, 
EW(H C ) = 7.6 ± 2.1 A, and EW(CaII H + H E ) = 16.7 ± 
2.2 A. 

Figure [4] shows D(CaII) measured in synthetic spec- 
tra for various metallicities as a function of ages of stellar 
populations. The synthetic spect ra were obtained based on 
the instantaneous-burst models (jBruzual fc Charlotl lioolh 
for metallicities Z = 0.005, 0.02, 0.2, 0.4, 1, and 2.5 Zq. 
The horizontal line shows the constraint from the observed 
D(CaII) at 1, 2, & 3 a upper limits. This diagram indicates 
that the inferred age of SST J1604+4304 depends on the 
metallicity; the 3 a constraint is <1,000 Myr for Z = 0.2Zq, 
<750 Myr for 0.4 z7 ,<6OO Myr for 1.0 Z ,and <400 Myr 
for 2.5 Z Q . 

The [OII]A3727 line flux is 5.5 ± 0.8 x 10" 17 ergs cm" 2 
s _1 . Its implications for the star formation rate (SFR) will 
be discussed later. 
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3.2 Analysis of the UV to near-infrared SED 
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Figure 5. Ax v versus metallicity, where Ax^ 



xl ~ xl(min). 



2.5Zq with the foreground Calzetti's dust screen. The solid line 
shows Ax^ for foreground dust screens with the Calzetti's extinc- 
tion law, the dotted line for foreground dust screens with the MW 
extinction law, the dashed line for internal dust models with the 
Calzetti's extinction law, and the dash dotted line for internal dust 
models with the MW extinction law. Note that \i was minimized 
for a given metallicity by changing ages f and E(B — V). 
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Figure 6. Ax^ maps for foreground dust screens (top) with 
2.5Zq instantaneous-burst models and internal dust models (bot- 
tom). The Calzetti's extinction law was used. Notice that internal 
dust models require much longer ages of stellar populations than 
those of uniform dust screens. 



We explore the nature of stellar populations and interstel- 
lar material in SST J1604+4304, by fitting model spectra 
of stellar population models to the observed SED from the 
B to 5.8um band. The 8.0um band is not included in this 
analysis, because emission from hot dust would be signifi- 
cant in this band. The ev olutionary synthesis codes, BC03 
l|Bruzual k, Cha riot 2003) are used to generate synthetic 
spectra of evolutionary stell ar popula t ion m odels. The ini- 
tial mass function (IMF) bv lChabrierl l|2003h (his table 1) is 
adopted with a lower and upper mass cutoff at 0.1 and 100 
Mq, respectively. 

We use two extinctio n laws , namely, the starburst ex- 
tinction curve bv lCalzettil (1200 ll) a nd the Milky way (MW) 



extinction curve by iDraine (2003). The starburst extinc- 



tion curve was derived from the data in the UV and op- 
tical wav elength range of local starbursts and blue com pact 
galaxies (|Calzetti. Kinney. fc Storchi-Bergm ann 1994). We 
consider two extreme cases of geometrical distribution of 
dust, namely, a foreground dust screen and internal dust 
model. In case of the foreground dust screen, dust is uni- 
formly distributed in the screen which is physically distant 
from star clusters, so we observe I°e~ T where r is the opti- 
cal depth of the dust screen and 7° is unextincted radiation. 
In the internal dust model, dust and stars are uniformly 
mixed and the model schematizes the situation in which 
dust is purely internal to the staburst site, so we observe 
J°(l-e- T(1 - 7) )/[r(l-7)], where 7 is the albedo. Note that 
we assume the isotropic scattering. 

ICalzettil (|200lh gives the extinction curve in a narrow 
range from 0.12 - 2.2 urn. To analyze multi-wavelemgth data 
from the UV to infrared, the extinction curve should be ex- 
tended t o wavelengths > 2.2 um. The analytical expression 
given bv lCalzettil j|200ll ) cannot be used for this extension, 
because the expression makes extinction values sharply drop 
beyond 2.2 urn . The Calzett i's extinction curve is almost 
identical to the lDraind l|2003l ) MW curve from 0.6 -1.6 urn. 
We thus extend the Calzetti's curve to longer wavelenths by 
adopting the MW extinc tion curve for w avelengths ^ 1.6 
urn. In addition, because ICal zctti (200 J) does not pro vide 
the albedo, we adopt the albedo given in iDraind {2003) for 
the Calzetti's extinction law with removing the dip in the 
albedo which is attributed to the familiar 2175 A feature. 

We start with instantaneous-burst models where all 
stars formed at once. The best model fit was searched for 
the observed SED by simultaneously optimizing the parame- 
ters, i.e., flux-scales, ages of stellar populations t, extinction 
E(B - V) intrinsic to SST J1604+4304, and metallicities of 
the stellar populations(Z). The ranges of the parameters are 
t = 0.01 - 6 Gyr, E(B — V) = 0.0 - 2.0 for the foreground 
dust screen and E(B — V) = 0.0 - 12.0 for the internal dust 
model, and Z = 0.005 - 2.5 Z®. 

The results are summarized in Figures [S] and [5] The 
best-fit was obtained for t = 40 Myr, E(B-V) = 0.8, and 
Z = 2.5Zp) with th e foreground dust screen combined with 
the ICalzettil l|200ll ) extinction law, resulting in a reduced 
chi square \t of 5.2. This xl is the minimum among all the 
models analyzed in this work. We thus refer to this minimum 
as xt( m i n ) = 5.2, and the increment Xv re lative to x^ jmin) 
is represented as Axt — xt ~ X^( m * n ) (|Avni| [l976V The 
approximate values of Axt are 1.15, 1.83, and 2.77 for 68%, 
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Figure 7. Instantaneous-burst models with an age 40 Myr for 
various metallicities. From the top to the bottom, the synthetic 
spectra for Z = 0.2, 0.4, 1.0, &2.5Zq are given. 



95%, and 99.6% confidence in this case(i.e., the degree of 
freedom is v — 10). 

The foreground dust screen with the MW extinction 
law is ruled out, because of the too large Axt(> 3). The 
internal dust models are ruled out either for two reasons: 
(l)inferred ages of the stellar populations are too large to 
be fit within the range allowed by the D(CaII) index; (2) 
Axv is too large, i.e., A%„ > 3 for the Z ^ Zq models and 
Axt = 1.8 for Z = 2.5Zq. For example, the best model for 
the internal dust geometry is obtained for Z = 2.5^0. As 
can be seen in the bottom panel of Figure [6j the inferred 
age is t — 0.7 - 0.8 Gyr. This age is not allowed because 
the age range must be < 400 Myr, corresponding to the 3<r 
D(CaII) limit for Z = 2.5Zq , as indicated in Figure[4] Figure 
[5] shows that models with a rich metallicity give better fit to 
the data than those with a low metallicity. The sensitivity 
to metallicity arises from the broad concave feature around 
rest 7000A, as shown in Figure [7] This feature is very weak 
in low metallicity spectra and becomes outstanding toward 
the high metallicity. 

We also examined models with the exponentially 
declining star formation rate, which is represented by 
T~ 1 exp(—t/T e ), where r e is the e- folding timescale and t is 
the time after the onset of the first star formation. r e ranges 
from 0.01 - 5.0 Gyr in a step of 0.5 dex and t from 0.01 - 
7 Gyr. The Calzetti's extinction law was used in the uni- 
form dust screen geometry. All models with Z Zq are 
rejected because they have Axt > 2. The D(CaII) index al- 
lows 2.5 Zq models to have t ^ 500 Myr, corresponding to 
3cr D(CaII) limit for exponentially declining star formation^. 
When t is fixed and the other parameters are changed, the 
smallest values of xt are obtained for E(B — V) — 0.8. In 
addition, t correlates well with r. Using this correlation, we 
found r ^ 1000 Myr. Applying another constraint Axt < 1, 
corresponding to the 1 a deviation, we find that r = 90 Myr 



4 For a given age of galaxies, D(CaII) is smaller in exponen- 
tially declining star formation than in instantaneous-burst, be- 
cause in exponentially declining star formation stars are continu- 
ously forming. 



and t = 200 Myr are allowed for the exponentially declining 
star formation models. 

The results are summarized in Table [2] In Figure [8] we 
plot the best-fit(i.e., A\ 2 = 0) instantaneous-burst model 
with Z = 2.5Zq and an age of 40 Myr. 



4 DUST EMISSION AND ENERGETICS 
4.1 Infrared luminosity 

As shown in Figure [8] dust emission is prominent at the 
24, 70, and 160u.ni bands. Because no data are avail- 
able longward of 160 u.m we fit the Arp 220 template 
to these data points to estimate the bolom etric infrared 
lumin os ity. The template is made based on [Carico et al.l 
(1990). Rigopoulou. Lawrence.fc Rowan-Robinson ~ifl996). 
and Klaas et al. 1 19971 ). In Figure 151 the upper dashed line 
represents the template which is fitted to the three data 
point at 24, 70 and 160 u.m, while the lower dashed line 
shows the template fitted to the two data points at 70 
and 160 \im. The former seems to overestimate the lumi- 
nosity, while the latter could underestimate the luminos- 
ity. We adopt the average of these two values for the real 
luminosity. Then, the 8-1000u.m bolometric luminosity is 
L ir = 1.78 ± 0.63 x 10 12 L Q , classifying SST J1604+4304 
as an ULIRG. Note that Li r includes the correction for the 
flux magnification by cluster CI 1604+4304, which is esti- 
mated to be 1.17. 



4.2 Energetics 

The stellar luminosity, which is inferred from the best 
instantaneous-burst model, is 8.9 x 10 11 Lq, accounting for 
only 1/2 of Li r . If star formation activity dominates the 
bolometric output, the significant part of the starburst site 
is hidden by very opaque dust extinction which we cannot 
observe at UV to near-infrared wavelengths, implying the 
dust distribution in the foreground screen is clumpy. Other- 
wise, 1/2 of Li r is fueled by AGN emission. 

The XMM - Newton X-ray 0.2-10 keV and VLA 
radio 20cm data listed in Table [T] do not support the 
AGN emission in SST J1604+4304. However, it should 
be noticed that the X-ray and radio observations are 
not deep enough to rule out a significant contribution 
from the AGN to the bolometric luminosity Lj r . Mid- 
infrared data are widely used to search for obscured 
AGNs, because these wavelengths are much less affected 
by extinction and can penetra te through the d usty atmo- 
spheres surrounding the AGNs (|Park et al]|2008l ). Using the 
f v (8u m) //„ (4.5/xm) versu s f v ( 2Aum) I fv(8um) AG N diag- 
nostic, llvison et ail (|2004T ) and I Pope et all |2008a|) showed 
that infrared galaxies having /„(8/im)//„(4.5/im) < 1.9 are 
starburst dominated in the mid-infrared. SST J1604+4304 
having f u (8fim)/ f v (4.5fim) — 1.22 ± 0.12 is thus starburst 
dominated. SST J1604+4304 has bluer mid-infrared col- 
ors th an those of ob s cured AGNs and t y pel Q SOs stud- 
ied by iPolletta et al] (|2008T ). IStern et al l (|2005l ) showed a 
[3.6/im] — [4.5/im] versus [5.8um] — [8.0/im] color-color dia- 
gram for identifying AGNs. SST J1604+4304 with [3.6/xm]- 
[A.bfxm]{Vega) = 0.35±0.08 and [5.8fim]-[8.Qfj,m](Vega) = 
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0.85 ± 0.12 lies on the edge of the region wh ich empir- 
ically separates AGNs from normal galaxies. lYan et al.l 
(2007) performed mid-infrared spectroscopy of infrared lu- 
minous galaxies at z ~ 2 with the sample selection crite- 
ria vfv{2&iim)/vf v {&iim) > 3.16. SST J1604+4304 hav- 
ing uf v (24:fim)/vf l , (8/im) = 1.6 does not meet the crite- 
ria of their sample. In conclusion, it is likely that mid- 
and far-infrared emission is starburst-dominated in SST 
J1604+4304, however we are unable to rule out the pres- 
ence of AGN emission. 

We compare the UV to far-infrared SED of SST 
J1604+4304 with those of three starburst-dominated galax- 
ies, namely, M82(local starburst), Arp220(local ULIRG), 
and MIPS J142824.0+352619(z = 1.3 HyLIRG) in Fig- 
ure |9] MIPS J142824.0+352619 with L ir = 3.2 x 
lO 13 Z/0 lacks any trace of AGN activity and has 
large PAH features and rich molecular gass es, indicat- 
ing starburst-dominated mid-infrared e mission (|Desai et all 
l2007llBorvs et all l2006lllono at al.ll2006l ). As seen in Figure 
[9] there is a striking similarity between SST J1604+4304 
and MIPS J142824.0+352619 in the overall SEDs from the 
UV to far-infrared. 



4.3 Dust mass 

The dust mass Mdust is related to the flux density f v as 
M dust /f v = D 2 /\K v Bu{T duat )] in the rest frame, where D 
is the distance to the source an d k v is the absorpti on cross 
section per unit dust mass(e.g., lOvabu et al. i 120091 ) . When 
a modified blackbody spectrum ' B u (Td U st) is applied, (3 
and Tdust should be used together as a set, because Mdust 
depends on K v B v (Tdust), but not on B v (Tdust) alone. To de- 
rive Mdust, we w ill use the du st model for the Milky Way in 
the far-infrared (|Drainell2003h in which k v oc v 2 or f3 = 2. 
Assuming (3 = 2, we obtain 32.5 K for the Arp 220 tem- 
plate and 35K for MIPS J142824.0+352619. This range of 
temperatures implies Mdust = 2.0 ± 1.0 x 10 s Mq for SST 
1604+4304, where the uncertainty is dominated by the error 
in the 160 urn flux. 
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Figure 8. UV to far-infrared SED for SST J1604+4304 along 
with the best-fit stellar population model with the Arp 220 
template. The solid line shows the unreddened stellar spectrum 
and the dotted line the stellar spectrum with the uniform dust 
screen E(B — V) = 0.8. This is the best-fit instantaneous- 
burst model for an age 40 Myr. The dashed lines represent 
the Arp 220 template fitted to three data points at 24, 70, 
& 160 |rm (upper) and two data points at 70 fc 160 yim 
(lower ) , where the template is made based on ICarico et al.l 
Il990h. iRigqpoulou. Lawrence. fc Rowan- Robinson! l| 19961 ). and 
Klaas ct al. I lll997l ). 



Table 2. Pure star formation with no AGN for SST J1604+4304 



Quantity 


Instantaneous™ 


E-decliningk 


E(B - V) c 


0.83 


0.79 


z(z Q r 


2.5 


2.5 


T e (Myr) e 




90 


t(Myr) f 


40 


200 


LI Ml 


28 


14 


Lir(L(7)) 


1.78 ± 0.63 x 10 12 


1.78 ±0.63 x 10 12 


M*(M ) 
Ax? 


6.4 ±2.3 x 10 10 


12.7 ±4.5 x 10 10 


0.0 


0.98 



5 DISCUSSION 

For simplicity, we temporarily assume that SST J1604±4304 
is purely powered by star formation and at the end we will 
discuss the impact if SST J1604±4304 is AGN-dominated. 



5.1 Stellar mass 

The stellar mass M* can be estimated as Li r / (L /M*) , where 
L/M* is the luminosity to mass ratio for the stellar popula- 
tion models. Thus, we have M, = 6 - 13 x 10 10 M Q at t = 
as listed in Table [2] If this stellar mass was constantly built 
up during 40 - 200 Myr, the star formation rate is SFR = 
600 - 1650 Mq yr _1 . Us ing the calibrati ons of SFR in terms 
of the [Oil] luminosity l|Kennicuttlll998l ). after dereddening 
the [Oil] flux, implies 13 - 500 Mq yr~\ Considering that 
we are only looking at 1/2 of the starburst activity in the 
the optical, we obtain SFR = 26 - 1000 Mq yr _1 , in good 
agreement with the previous estimate. 



a Instantaneous-burst models where stars form all at once. 

6 Exponentially declining star-formation models with SFR 

oc Te^e~ t l T ", where t is the time after the onset of the initial 

star f ormation and r e is the e-folding timescale. 

c The Calzctti ( 2001) extinction law is used in a foreground dust 

screen. 

d Metallicity of the stellar populations. 
e e-folding timescale for SFR. 

* Time after the onset of the initial star formation. 

9 Luminosity to mass ratio in solar units, where M* is a stellar 

mass at t = 0. 



5.2 Comparison with other infrared galaxy 
populations 

SST J1604±4304 is a dusty, super metal-rich, young galaxy 
at z = 1.135 with Eb-v = 0.8, metallicity 2.5 Zq, and 
an age 40 - 200 Myr. The age is comparable to young, 
optical-selected populations like Lyman-a emitters (LAEs) 
and Lyman break galaxies (LBGs). Studying the z ~ 3 
LAE sample, iGawiser et al.l (|2006l ) inferred ages of the or- 
der of 100 Myr with stellar masses 5 x 10 8 A/q and no 
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Figure 9. UV to far-infrared SED for SST J1604+4304 
compared with those of starburst-dominated galaxies, 
M 82(dot local starburst), Arp 220(solid local ULIRG), 
MIPS J142824.0+352619(das/i z = 1.3 HyLIRG ). M82 



data come from [ Willner et al. | (119771), iDale et all 
iHughes. Gear, fc Robs onl (119941). iLeeuw fc^obsorT 



2001) 



2009) 



and IWeedman*"*CIiarni andaris. &: Zeza si (1200411. Arp 220 from 
ICarico et al. I dl99Ch. Rigopoulou, Lawrence, & Rowan-Robinson 
dl996h, and iKlaas et al.l dl997T). and MIPS J 142824.0+352619 
from fpesai et a ll ( 12007^ 1 and lBorvs et al. (2006). 



dust extinction. iPirzkal et al.1 (|2007t > showed 4 < z < 5.7 
LAEs had ver y young ages a few Mys w ith stellar masses 
10 6 - 1O 8 M . ISawicki fc Yee et al.l (|l998l ), studying z > 2 
LBGs, found that their LBGs were dominated by young stel- 
lar populations with ages < 200 Myr with stellar masses 
several x10 9 Mq and moderate dust extinction typically 
E(B - V) ~ 0.3. The z = 2 - 3.5 s ample of LBGs by 
IPapovich, Dickinson, fc Ferguson! (2001) has ages of 30 Myr 
to 1 Gyr with stellar masses 10 9 - lO n M with E(B - V) 
= 0.0 - 0.4. SST J1604+4304, having a stellar mass 6 - 13 
x10 10 Mq, is more massive and more dusty than these op- 
tically selected galaxies. 

SST J1604+4304 having /„(24^m) //„(#) > 1000 
meets the criteria for 24(j.m-selected high-redshift (z ~ 2) 
dust-obscured galaxies (DOGs). DOGs and submillimeter- 
selected galaxies (SMGs) have extremely large Li r /Ls 
which cannot be accounted for by r edshifting l ocal U LIRGs 
dDev et alll2008l|Pope etalll2008bl ). iDev et all (|2008h found 
nearly all of DOGs are ULIRGs with L IR > 10 12 L Q and 
more than half of these have Lib. > 10 13 Lq. The infrared 
luminosit y of radio-identified SM Gs ranges Lib — 2 x 10 11 — 
1O 14 L (|Chapman et all 120051 ). A SF R 1000 M^yr~\ 
which corresponds to Lib = 5 x 10 12 Lq (|Kennicuttl fl998). 
builds up a stellar mass of 10 11 Mq in 100 Myr. Thus, it is 
considered that DOGs and SMGs are progenitors of present- 
day massive galaxies. If ages of stellar populations in DOGs 
and SMGs are derived, as we did for SST J1604+4304 by 
measuring the stellar absorption lines, the evolutionary link 
of DOGs and SMGs to normal ULIRGs and optically se- 
lected galaxies will be better understood. 



because low-mass stars had not evolved to form dust in 
the expanding envelope. This may be the reason why our 
best model results in poor goodness of fit xl(min) — 5.2. 
The Calzetti's extinction law may not be applicable to such 
young galaxies, because dust grains in local starbursts would 
be dominated by those formed in the outflows of low-mass 
stars and not by SNe II dust. In this regard, it is of great 
interest to use the SED of SST J1604+4304 to test extinc- 
tion laws such as those calculated by Hirashita et al. (2005, 
2008) that are pred icted for dust formed in SNe II ejecta 
jNozawa et al.ll2003t ). 

Our best model is the instantaneous-burst model with 
an age 40 Myr with a stellar mass 6.4 ± 2.3 x 1O 1O M 
at t = 0. The metallicity is 2.5 Zq and the dust mass is 
M dust = 2.0 ± 1.0 x lO 8 Af . Until 40 Myrs after the on- 
set of the starburst, stars with a mass > 8Mq died with a 
SN explosion. The number of stars evolved into SNe in 40 
Myr is 7.6 x 10 8 or 1 9 SNe per yea r . Ado pting the nucle- 
osynthesis models by iNomoto et al. (2006), the total yield 
is 0.09 - 0.11 and the metal yield is 0.016 - 0.020 for the 
Chabrier IMF0 The gas and stetllar metallicities can be 
estimated applying simply analytical models based on the 
instantaneous recycling approximation. The gas metallicity 
is expressed as Z g (t) = y(t)ln(l/ f(t)) where f(t) is the ratio 
of the gas mass to the galaxy mass, and y(t) is the ratio of 
the rate at which the metal is produced by the events of nu- 
cleosysnthesis and ejected into the interstellar gas to the rate 
at which hydro gen is removed from th e interstellar gas by 
star formation (Sc arle fc Sargendll972l ). The stellar metal- 
licity is derived from S(Z)/Si = (1- f? /Zx )/(l- fi) and 
Z, = J Q Zl ZdS{Z)/ J Zl dS(Z), where S{Z) is the total mass 
of stars born up to a time when the metal abundance reached 
to the valu e Z, and Si,Zi,fi deno te S(ti), Z(ti), /(ti), re- 
spectively dPagel fc Patch ctt 1975). In our case, t\ is a time 
when the stellar metallicity reached to Z* — 2.5Zq. When 
the gas mass is reduced to / = 0.3, we have Z, — 2A—2.7Zq 
and Z g — 4.5 — 5.1 Zq . The inferred galaxy mass is 1.2 x the 
stellar mass at t = 0, because the total mass ejected from 
stars at the events of SN explosions into the interstellar gas 
is 0.1 per unit stellar mass at t = 0. The total gas mass is 
2.3 x 1O 1O M . The gas-to-dust mass ratio is then 120 ± 73, 
which is comparable t o the value obta ined in a z = 1.3 
HyLIRG (198 ± 53) bv llono at all (|2006l ). 

The dust production rate is O.24±O.12M per SN. This 
is consistent with the numerical results of dust formation in 
SNe; for progenitor masses ranging 13 - 30 Mq, 2 - 5% of 
the pr ogenitor mass is loc ked into dust grains at SN explo- 
sions (|Nozawa et al.ll2003l ). 20 - 100% of which is destroyed 
by processing through the collisions with the reverse shocks 
resulting from the interaction of SN ejecta and with the am- 
bient medium (|Nozawa et al.l [20071 ). 



5.4 Foreground dust screen 

As demonstrated in Section 3.2, a foreground dust screen 
geometry is plausible for SST J1604+4304. This indicates 



5.3 SNe for dust and metal production 

In SST J1604+4304 with an age 40 - 200 Myr, most of 
the dust would be formed by type II supernovae (SNe II), 



5 The original yields are given for the Salpeter IMF. These are 
converted to the Chabrier IMF by simply multiplying 1.64 to 
adjust the difference in the mass fraction of stars from 8 - 100 
M & to those from 0.1 -100 Mq. 



10 K. Kawara et al. 



that dust is depleted in the starburst site and the stellar 
spectra are reddened by the foreground dust screen which 
enclosed the starburst site. Let's suppose a spherical shell 
with a radius I for the dust distribution, then we observe 
/„ = 4T v B v (Tdust)ir(l/D) 2 in the rest-frame, where the 
thickness of the shell is t v which is 0.0074 at re st 80[i.m corre - 
sponding to E(B-V) =0.83 with R v =4.05 (|Drainell2003l ). 
D is the distance from the observer to SST J1604+4304, 
and Tdust = 32.5 — 35K. Converting this relation into the 
observer's frame, we obtain the radius I — 4.5 - 5.5 kpc or 
0.56 - 0.65". This is comparable to the optical size of SST 
J1604+4304 which is approximately 1.2" x 0.5"in the J 8 i4 
band. This supports the view that the galaxy is surrounded 
by the shell. Further support is the fact that there are no 
significant color gradients across the galaxy in the HST im- 
ages. 

5.5 Evolutionary link to other galaxy populations 

As discussed by ICalzettil (|200ll ). starburst environments 
are rather inhospitable to dust; dust grains in the star- 
burst site can be transported to a larg e distance in a rel- 
atively short time by ra diation pressure l|Ferrara et al.lll99ll ; 
IVenkatesan et al1l2006l ). as well dust grains that formed at 
SN explosions are processed and evolve in SN remnants 
l|Nozawa et al.l l2007) - small size grains are quickly destroyed 
in SN remnants by sputtering. Thus, a cavity-shell struc- 
ture is a natural geometry for the dust distribution in star- 
forming galaxies; the starburst site is inside the cavity where 
dust is depleted, and the opaque dust shell is surrounding 
the cavity. The shell is observed as a foreground screen. Such 
foreground screens were found in local starburst galaxies 
llCalzetti. Kinney, fc Storchi-Bergmannlll994l ; iMeurer et al.l 
119951 ). In the foreground screen of SST J1604+4304. the dust 
would be unevenly distributed; we are only looking at part 
of the starburst site at UV and optical wavelengths through 
relatively transparent holes with E(B — V) = 0.8, and the 
other part is completely obscured at these wavelengths. This 
explains that Li r is two times greater than the stellar lumi- 

nosi ty derived from the broa d-ba nd SEP analys is. 

IVenkatesan et al.l (120061 ) and iNozawa et~atl (|2007l ) sug- 
gest that dust created in the first SN explosions can be 
driven through the interior of the SN remnants and ac- 
cumulated in the SN shells, where second-generation stars 
may form in compressed cooling gases. Hence, galaxies may 
be observed as dust-free objects at the very beginning. 
When much dust formed in SN ejecta, the starburst site 
would become completely obscured by dust at UV to near- 
infrared wavelengths. This stage would correspond to ex- 
tremely dusty objects like DO Gs and SMGs. Then, the 
galact ic winds would turned on. iHeckman. Armus. fc Milevl 
l|l990l ) discussed that the galactic winds or superwinds fre- 
quently observed in starbursts and ULIRGs will sweep out 
any diffuse interstellar matter from the starburst site. Such 
galactic winds will turned on when the thermal energy of 
the gas heated by SN explosions exceeds the gravitation- 
ally binding energy of the gas. According to the models for 
chemi cal evolution of elliptical galaxies bv lArimoto fc Yoshiil 
galactic winds turn on later in more massive galax- 
ies. The onset of a galactic wind is 350 Myr for a 10 1 1 Mq 
galaxy, and the metallicity increases up to more than the so- 
lar value. These predictions seem to agree well with the age 



and the metallicity observed in SST J1604+4304, in which 
the galactic wind would just turn on and dust grains in the 
opaque shell are pushed and moved outward, creating par- 
tially transparent holes in the shell. This picture is consis- 
tent with the fact that dusty objects are deemed to be more 
abundant in massive galaxies than in less massive galaxies, 
although mid- and far-infrared observations are strongly bi- 
ased to IR-luminous objects. This observational trend is ex- 
pected, because galactic winds turn on later and thus the 
dust shell becomes more opaque for more massive galaxies, 
resulting in longer time of the dust-obscured stage. This may 
be the reason why DOGs and SMGs are so infrared-luminous 
and massive. 

5.6 Impact by an obscured AGN 

What would happen if the missing 1/2 of the bolometric 
luminosity comes from an AGN? In this case, the stellar 
mass ranging 0.1-8 Mq is reduced to one half of that 
for the pure star formation. The number of SNe is also re- 
duced by the same amount, while the dust mass does not 
change. Thus, the dust production rate is increased by a fac- 
tor of 2, i.e., 0.48 ± O.24M per SN. T his is within a range 
consistent with the model predictions dNozawa et ail 120031 ; 
INozawa et al.ll2007l ). The metallicity is the same as derived 
for pure star formation, because the SN ejecta mass is also 
reduced by a factor of 2. If the AGN shines at the Eddington 
rate, the mass of the AGN is 2.5 x 1O 7 M . 
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